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1.0  Summary 


The  main  effort  this  year  is  to  understand  the  physics  of  the  hybrid  plume  and  heating 
of  the  plasma  in  a  tandem  mirror  device.  In  the  theoretical  area  three  tasks  have  been 
accomplished:  (1)  .A  major  advance  has  been  made  in  the  development  of  a  numerical 
method  for  solving  the  time  dependent  fluid  equations  for  both  ions  and  electrons  in  the 
plasma.  Realistic  physical  properties  and  nonuniform  magnetic  fields  were  considere<I.  (2) 
-A  mathematical  model  has  been  formulated  for  studying  the  radio  frequency  heating  of 
the  hot  plasma  in  cylindrical  geometry.  (3)  .An  analysis  on  the  performance  characteristics 
of  the  hybrid  plume  rocket  has  been  carried  out. 

Under  the  funding  of  the  instrumentation  program  some  effort  has  been  dedicated  to 
the  construction  of  a  tandem  mirror  device.  The  component  fabrication  began  in  May  after 
several  months  of  contract  negotiations  since  authorization  of  the  fund,  and  is  progressing 
well. 

To  obtain  a  quick  but  crucial  understanding  of  the  hybrid  plume  concept,  a  steady 
state  fluid  model  code  for  studying  the  interaction  of  hot  plasma  with  neutral  gas  jet 
was  developed  in  past  years.  In  that  code  a  uniform  magnetic  field  was  assumed.  It  was 
found  that  the  transition  region  has  been  formed  and  an  exhaust  duct  can  be  protected 
by  surrounding  the  hc)t  plasma  with  a  supersonic  gas  jet.  Therefore  the  simplified  version 
of  the  computation  code  has  fulfilled  the  purpose  of  verifying  the  hybrid  plume  concept. 
To  further  understand  the  detailed  flow  dynamics  of  the  exhaust,  a  time  dependent  fluid 
model  for  plasma  in  a  general  magnetic  field  configuration  has  to  be  solved  with  realistic 
physical  assumption.  To  develop  such  a  model  was  the  task  undertaken  this  year.  .After 
overcoming  many  difficult  numerical  instability  problems,  a  fundamental  two  dimensional 
time  dependent  code  is  now  operational.  More  importantly,  the  plasma  is  now  treated  as 
two  fluids,  the  ions  and  electrons,  instead  of  one  fluid  used  in  the  previous  code.  However, 
many  difficult  problems  have  to  be  worked  out  in  order  to  use  the  code  in  a  parametric 
study.  The  outflow  boundary  has  to  be  properly  treated.  The  neutral  gas  jet,  atomic 
reactions  and  radiation  have  to  be  included.  Solving  the  outflow  boundary  problem  will 


be  the  next  major  task  for  the  coming  year. 


Previously  a  slab  model  for  studying  the  ion  cyclotron  range  of  frequence  heating 
was  developed.  It  was  shown  that  the  heating  efficiency  can  be  impnaed  by  optirni/ing 
the  antenna  loading  impedance.  The  loading  impedance  on  the  antenna  was  found  para¬ 
metrically  dependent  on  the  magnetic  field  strength,  the  distance  between  the  antenna 
and  vessel  wall,  the  plasma  density  and  the  dimension  of  the  antenna.  Currently,  a  more 
accurate  and  realistic  cylindrical  mixlel  for  the  power  absorption  in  the  plasma  has  been 
fi-irmulated  analytically.  In  order  to  learn  the  heating  mechanism  and  to  optimize  the  effi¬ 
ciency,  we  would  like  to  obtain  an  explicit  parametric  dependence  of  the  power  absorption 
in  the  plasma  on  machine  and  plasma  parameters. 

For  the  construction  of  the  tandem  mirror  device  the  central  cell  coils  and  one  mirror 
coil  have  been  delivered.  The  liquid  nitrogen  coil  casings  are  in  the  process  of  being 
fabricated.  The  fabrications  of  the  vacuum  vessel  is  progressing  well.  The  assembly  of  the 
device  is  expected  in  .November. 


2.0  Concept 

It  woulil  be  useful  to  reintroduce  the  concept  in  this  section  and  to  investigate  the 
physical  problems  in  the  next  section.  .A  hybrid  plume  rocket  is  one  where  the  macroscopic 
properties  of  the  exhaust  fluid  (i.e..  temperature,  densit}'.  fluid  velocity,  etc.)  exhibit  a 
drastic  variation  from  one  radial  position  to  another  1,5  .  In  a  hybrid  plume  rocket,  the 
temperature  at  the  core  of  the  exhaust  is  extremely  high;  whereas,  near  the  edge,  it  is  quite 
low.  In  fact,  the  \ariation  in  temperature  from  the  plume  centerline  to  the  outside  can  be 
many  order  of  magnitude,  while  the  exhaust  fluid  at  the  center  is  a  fully  ionized  plasma  (at 
.5  to  1  ke\’).  at  the  edges  it  actually  becomes  a  neutral  gas  (at  .1  to  .3  e\').  .\s  a  result,  the 
cool  flow  at  the  edges  relaxes  the  heating  constraints  on  the  nozzle  material  while  the  hot 
central  core  preserves  and  even  enhances  the  allowable  power  densit\’  and  specific  impulse 
of  the  device.  The  concept  has  a  potential  application  in  rocket  propulsion. 
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Figure  1  shows,  the  tandem  mirror  device  with  a  hot  [)lasma  column.  One  could 
envision  a  very  hot  plasma  exhaust  surrounded  by  an  annular  layer  of  insulating  gas. 
The  hot  plasma  is  generated  by  a  two  stage  process  involving  an  initial  source  of  cool 
dense  plasma  injected  axially  into  a  straight  tandem  mirror  device  from  the  left  end.  Once 
magnetically  confined,  the  plasma  can  be  heated  to  the  ke\'  regime  by  microwave  injection. 
The  plasma  then  escapes  through  one  end  into  an  exhaust  nozzle  on  the  right  end.  The 
hybrid  plume  is  formed  slightly  downstream  of  the  nozzle  throat,  as  the  hypersonic  annular 
gas  layer  is  injected  and  a  stable  bountlary  is  established. 

The  rocket  propulsive  efficiency  q  is  in  free  space  and  in  the  absence  of  gravitational 
effects  is  qp  ---  I'y^V  in  terms  of  vehicle  velocity  v  and  exhaust  velocity  u.  Optimum 
efficiency  occurs  when  r  -  u.  The  hybrid  plume  would  produce  a  variable  specific  impulse 
I,p  which  allows  the  engine  to  be  optimally  tuned  for  all  phases  of  flight.  The  terminal 
velocity  of  the  hybrid  is  nearly  ten  times  the  initial  velocity  when  reaching  burnup. 

3.0  Physical  Problems 

The  physics  (,)f  high  temperature  hvdmgenic  plasmas,  interacting  with  hypersonic 
gas  jets  are  not  well  understood.  .Analytical  and  inimerical  modeling  of  these  interac¬ 
tions  requires  the  simultaneous,  self-consistent  solution  of  particle,  momentum  and  energy 
transport  equations  for  each  particle  sp<  cies.  Therefore,  the  crucial  area  of  investigation 
pertaining  to  the  phy.sics  of  hybrid  plumes  is  the  numerical  modeling  of  the  plasma.gas 
interaction  in  the  exhaust. 

The  numerical  modeling  results  have  to  be  tested  experimentally.  The  plasma  con¬ 
ditions  achieved  in  the  fiasion  program  have  exceeded  the  need  fur  propulsion  purposes. 
However,  the  devices  used  for  fusion  experiments  are  usually  large  because  the  purpose  is 
aimed  at  power  reactors.  To  find  a  fight  and  compact  dev  ice  suitable  for  prr<pulsion  has 
to  be  carefully  investigated.  The  basic  plasma  physics  for  the  compact  device  and  control 
of  plasma  flow  are  not  known  and  need  extensive  experimental  study. 
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4.0  Time  Dependent  Fluid  Modeling 


As  discussed  in  the  suinniiiry.  the  parametric  study  obtained  by  the  code  developed 
previously  has  shciwn  the  viability  of  the  hybrid  plume  rocket  .  1  he  interactions 

between  the  neutral  gas  jet  and  plasma  have  been  demonstrated  to  form  the  transition 
region  requiretl  to  form  the  hybrid  plume.  .Also,  witli  \arying  the  amount  of  neutral  gas 
introduced  into  the  s\  stem.  one  can  modulate  the  central  core  temperature.  I  his  is  an 
mifiortant  parameter  in  controlling  the  plasma  flow.  1  he  next  critical  issue  to  address  is 
the  detailed  flow  dynamics  of  the  hybrid  plume  concept.  Specifically,  whether  the  plasma 
flow  will  deiach  from  the  nozzle  and  proiluce  usable  thrust  at  high  I^p.  Preliminary  results 
show  that  the  neutral  gas  jet  can  provide  the  momentum  needed  to  cause  the  plasma  to 
form  closed  magnetic  flux  surfaces  which  will  flow  out  the  exhaust  duct  like  a  ring  and 
eventually  disconnect  from  the  main  field  line.  This  phenomena  will  be  referred  to  as  field 
line  disconnection  in  this  write-up. 

The  previous  code  used  a  cylindrically  symmetric  duct  with  the  plasma  flowing  axially 
and  diffusing  radially.  I’his  assumption  was  used  to  understand  the  principle  mechanisms 
involved  with  the  hybrid  plume  concept  and  to  study  the  feasibility  of  the  hybrid  plume 
concept.  The  code  was  not  designed  to  specifically  stiidy  the  disconnection  process.  The 
n.agnetic  field  was  assumed  constant  in  titru'  and  essentially  constant  in  space  (a  modified 
diffusion  coefficient  was  used  to  model  the  weakly  di\erging  magnetic  field).  The  code 
calculated  density,  n.  temperature  T,  and  axial  vehnity  \  I  he  radial  velocity  was  not 
needed  because  of  the  diffusion  model  used.  This  limited  the  interactions  fx-tween  the 
neutral  gas  and  plasma  to  the  diffusion  time  scale.  In  reality,  the  plasma  will  physically 
flow  radially  along  the  diverging  field  line.  In  a  real  nozzle,  1  j  will  be  important.  Depending 
on  the  strength  of  the  interactions  between  the  neutral  gas  jet  and  the  [)lasma.  the  [dasrna 
may  cling  to  the  field  lines  and  actually  reverse  directions  as  the  magnetic  field  bends 
around  the  edge  of  the  nozzle.  .\lso.  the  <lynamics  of  the  flow  ran  alter  the  local  magnetic 
field,  scj  the  field  would  be  non-stationary.  ffius  to  acccirately  model  the  h\brid  plume, 
and  the  separation  of  the  plasma  from  the  field  lines,  one  must  improve  the  numerical 
model  significantly. 


A  simple  mass  and  momentnm  balance  argument  shows  that  the  neutral  gas  jet  ran 
provide  the  momentum  needed  to  drive  the  disconnection  and  separate  the  plasma  flow 
from  the  field  lines.  Since  the  neiitral  gas  jet  is  located  near  the  wall  of  the  exhaust,  the 
interacti<'n  region  will  form  df)wnstream  from  tfiis  point.  I  he  magnetic  field  will  be  falling 
off  rapidly  in  this  region,  fhus.  the  momentum  required  in  the  neutral  gas  jet  can  be  small 
while  still  tlri\ing  the  reconnectit)n  pr<»cess.  .A  potential  sequence  is  shown  in  fig.  2.  The 
neutral  gas  Jet  flows  axially  down  a  duct  immersed  in  a  di\erging  magnetic  field.  .As  the 
neutral  gas  jet  becomes  ionized,  its  resistivity  and  flow  profile  modify  the  magnetic  field 
and  the  remaining  plasma  flow.  This  allows  the  flow  to  disconnect  from  the  nozzle  and 
continue  down  the  duct.  I  his  detachment  mechanism  is  presently  under  study  and  will  be 
one  of  the  new  modifications  to  the  present  code. 

Our  current  effort  has  been  to  improve  the  computational  model  to  include  a  more 
realistic  magnetic  field  geometry.  Further  refinements  in  a  physics  of  the  interactions 
between  the  neutral  jet  and  plasma  will  also  be  included.  In  an  actual  machine  it  is  not 
obvious  that  a  true  steady  state  will  be  obtained  (periodic  operation  can  be  envisioned). 
The  new  model  will  be  time  dependent,  two  dimensional  in  space,  and  three-dimensional 
in  velocity  and  electric  and  magnetic  fields,  for  fully  ionized  plasma  the  i(.)n  and  electron 
tlensity  are  approximately  equal.  These  two  species  of  particles  have  different  mobility 
and  -hould  be  treated  as  separate  fluids  in  an  open  system  like  our  case.  The  complexity 
i.'  double  that  of  the  single  fluid  treatment.  However,  the  electron  is  much  lighter  than 
the  ion.  It  is  reasonable  to  make  the  zero-electron  mass  assumption.  .A  zero-electron-mass 
h\ brid  code  is  being  developed.  It  is  a  two-fliiid.  two  dimensional  (r  and  c).  quasi-neutral 
code.  This  code  can  adequately  model  the  complex  flow  and  evolution  on  the  hybrid  plume 
dynamics. 

With  zero  mass  the  electron  momentum  equation  reduces  to 


VpTp  lie  X  E 
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w  here  E  is  the  electric  fu-ld  \  ector.  p  is  the  particle  density.  T,  is  the  electron  temperature,  c 
is  th  '  charge  of  an  electron,  ii,.  is  the  electron  velocity  vector.  E  is  the  magnetic  field  vector. 


q  is  the  resistivity  tensor,  and  J  is  the  current  density  vector.  The  fluid  equations  for  ions 
are  the  same  as  reported  previously,  (  sing  the  Darwin  limit  of  .Maxwell's  equations  and 
also  the  magnetic  vector  potential  A.  the  electron  momentum  equation  can  be  rewritten 
in  a  form  for  the  time  evolution  of  A, 

A  =  c— ^  -UeX(VxA)-— q-VxVxA.  (2) 

ep 

One  actually  only  needs  the  theta  component  of  the  above  equation  and  one  can  obtain 
both  Br  and  5..  .A  similar  transformation  results  in  an  equation  for  the  time  evolution  of 


The  electron  fluid  velocities  are  calculated  according  to  the  currents.  In  other  words. 
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,\11  that  is  needed  to  complete  the  picture  are  the  fluid  quantities,  and  they  are 
calculated  using  the  normal  fluid  equations  for  mass,  velocity,  and  temperature. 

The  gcial  of  the  current  development  will  be  the  abilit\  to  model  the  detailed  behavior 
of  the  plasma/gas  flow  in  a  complex  magnetic  field.  This  is  important  because  the  jilasma 
will  tend  to  cling  to  the  field  lines  (which  diverge),  while  the  neutral  gas  will  tend  to  go 
straight.  The  bulk  flow  will  be  some  complicated  balance  between  the  two.  I  he  higher 
the  plasma  temperature,  and  thus  higher  potential  /,p.  the  more  the  plasma  flow  will  be 
frozen  the  field  lines.  To  obtain  realistic  thrust  predictions,  detailed  characteristics  of  the 
flow  field  must  be  determined.  The  code  solves  for  the  e\olution  of  mass,  momentum, 
'■nerg\',  and  magnetic  field.  The  plasma  equations  are  solved  explicitly,  with  the  coupled 
field  equaticjns  solved  implicitly. 


5.0  Plasma  Heating 

Progress  made  in  ICRF  heating  this  year  includes  the  theoretical  study,  the  r 
absorption  of  the  plasma,  and  design  of  the  power  systems  for  ICRF  heating.  .\n  entii 
neering  design  of  the  rf  power  transmission  systems  was  performed  with  the  assistance  of 
the  radio  frequency  power  transmitter  expert  from  the  .Alcator  group. 

On  the  theoretical  side,  further  work  was  done  in  the  optimization  of  antenna  heatiim. 
Last  year,  power  absorption  calculations  were  performed  in  a  slab  Cartesian  geometrx. 
While  the  results  obtained  using  Cartesian  geometry  proved  insightful,  they  are  rmt  \ery 
physically  applicable.  The  geometry  of  the  device  is  such  that  the  free  space  wa\ (‘length 
of  the  RF  waves  generated  A  are  much  larger  than  the  physical  size  of  the  plasma  and  or 
the  device  a.  Since  the  Cartesian  approach  means  that  the  plasma  is  in  effect  of  infinite 
size  such  that  the  geometric  effects  of  its  cylindrical  shape  do  not  affect  the  results,  it 
means  that  we  need  A  a.  which  is  not  the  case  here.  So  while  the  calculations  obtained 
using  slab  geometry  are  useful  in  that  they  yield  intuitive  results,  they  are  not  a  very 
effective  model  of  the  physical  conditions  of  the  e,xperimenl.  This  year,  calculations  using 
the  cylindrical  model  are  performed  and  are  being  examined.  Collisions  are  included  to 
examine  what  elfects  collisions  may  have  on  the  heating  of  the  plasma.  In  the  future,  finite 
temperature  effects  will  also  be  introduced  to  better  model  a  laboratory  plasma. 

The  electromagnetic  field  structure  inside  a  cylindrical  plasma  was  calculated  and 
examined  to  maximize  the  power  injected  into  the  plasma  by  the  antenna. 
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I  rorii  Maxwell  s  equations,  the  fields  inside  a  cylinder  can  be  calculated  to  be: 


Region  1.  (Plasma  Region) 


//j  (  r  )  // 1  <I nlr  *'1^2*^n2r 


^-'2*^ri2r  *  y7-^l^nl 


Region  II.  i  \  acuiim  I^egion) 


//^(r)  ^  //„  Inr  '  Kjit  ■ 


E-(r)  =  E  /„^  ^  E  Kr,r  , 
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The  coefficients  of  the  fields  can  be  found  subject  to  the  following  boundary  conditions: 


(1)  .At  r^c. 
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(2)  .At  r-b 


Eeij  =  Eeijj  , 


Ezii  E-iii 


Hzii  =  H-iii  -  Jso 


flem  ^  Hen  -  J, 


f3)  .At  r  -a 


After  chlculations  were  perfurnied  to  find  the  field  coefficients,  the  power  injectetl  int' 
the  plasma  can  be  evaluated  by 
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Define  elements  of  the  dielectric  tensor  (including  collisions  of  collision  frequency  Uc): 
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^0  3  , 


-Xc’^^bci^  -  1)  _ 

,  ,  /  ^  ^±1  !  A 

-a  c  ^  ri2a  O/  ♦-'n2a*-'Qc' 


/,  ,  A' 

^  _  A  I-/  / 1  •  -QC  ' '  na  , 

•^2  =  iA(,rA„„(l  -  — - )  , 


)C  -  Q5fc^(n^  -  5) 

k^'PsnKr.c^^b'b 


The  results  obtained  from  the  analysis  yields: 
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Ee{r  -  a) 


'I'o-'to 


■^*61(^0 
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nAb'bEr.c'P^E 
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CJnla^t  '^b'h  KncV^ 


j  ,nk.Qo  k.aSk^,{2S  ~  n'^)A2  s 

akg  CJ„ian/\yt,f<nc'P'j 


The  radial  power  flux  is  found  by 


E  X  H* 


=  Esh: 


Assuming  a  loop  antenna,  the  formula  for  the  power  flux  through  the  plasma's  edge 


I'r  — a)  is  found  to  be: 


f^'O'^nSa*  ■r  2  / 


kz'HDn 


+  aPCi^  ^  D^ktm  -  klS  -  %kl,) 


-EuA\ 


N  , _ kz^g^S-  n])^  D 

'if  ^  ^  VlJn2aJn2a’  LJn2a^,.bKncV3P{^-^  ~  D-^kl){k]  -  klS  ~ 


kl^S[2S~nl)  \ 

^  2(<--^^2Jn2a  k^.'l'n Dokl  S  ^ {2S  -  n])[2S'  -  nl] 


.'i>Jn2aJn2.  Jr.2aPH'i>^  ^  D^-k^^}{ki  A'o  5  -  J  n£*  J.2, 


The  real  power  flux  absorbed  by  the  plasma  can  then  be  found  by  taking  the  real  part 
of  the  radial  Poynting  flux,  Re(Sr).  The  total  power  into  the  plasma  can  then  be  found 
by  integrating  the  Poynting  Flux  with  respect  to  k^. 


In  future  work,  the  cylindrical  analysis  will  be  extended  to  warm  plasmas  so  that 
a  more  realistic  model  of  a  laboratory  plasma  can  be  obtained  and  compared  to  future 


experimental  results. 
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G.O  Pcrfuriiianc'c  Characteristics  of  tiu'  Hybrid  Pliiinc  Rock(?t 


This  discussion  focuses  on  tiie  basic  physical  la\^s  governing  tlie  uperatiuii  and  <i[)ti 
niization  of  the  hybrid  plume.  It  assumes  a  suitable  power  smirce  can  be  pro\  ided  wdiicli 
can  deliver  a  power  kwel.  P,  to  the  plasma.  .Microwa\e  and  other  rf  efficiencies  are  not 
discussed  here  as  they  are  experimental  cpiantities  which  need  to  fie  measured  In  a  real 
ilevice.  .\s  such,  these  \alues  neitlier  add  nor  subtract  to  tfie  ph\sics  to  be  discussed  here. 
W  itfi  these  considerations,  the  f>asic  momentum  and  energy  relatioiisfiip'  as  thev  would 
apply  to  the  hybrid  plume  engine  will  not  f)e  explored. 

6.1  Basic  Assumptions 


The  basic  momentum  balance  for  a  rocket  in  free  space  and  assuming  no  gravitational 
effects  can  be  e.. pressed  as 

,  um' 

v'  --  ^  -  (-IT) 

m 

where  v'  and  m'  are  time  derivations  of  the  vehicle  velocity  with  respt-ct  to  a  fixed  point 
(i.e.,  the  Earth)  and  rocket  total  mass  respectively,  u  is  the  exhaust  velocity  as  measured 
in  the  frame  of  the  rocket. 

The  rocket  exhaust  power  flow.  P,  can  also  be  expressed  as 

O 

,  U" 

P  -  m  -  .  I  IS) 

.\lso.  the  various  masses  can  be  related  a  follows: 

m{t)  =  morrifit)  (19) 


where  rriQ  is  the  total  mass  of  rocket  and  fuel  at  t  -  0.  and  mf(j)  i.s  the  total  mass  of  fuel 
exhausted  tlirough  the  nozzle  at  time  t. 


We  not  define  the  propulsive  efficiency  of  the  engine,  qp.  as  the  fraction  of  the  total 
power  available  which  is  converted  to  vehicle  kinetic  power.  That  is. 


np 


(oO) 


I.d 


where  T  is  the  instantaneous  thrust,  m  u  and  the  quantity  m'(v  -  u)~  '2  represents  the 
residual  power  which  is  carried  out  by  the  exhaust.  In  terms  of  the  velocity  ratio  (r/u), 
Fiq.  (50)  can  be  recast  as  follows: 


2 


l-(  =  ) 


lol) 


liquation  (51)  is  shown  plotted  in  Fig.  5  as  a  function  of  the  velt)city  ratio.  It  can  be  seen 
that  the  optimum  occurs  when  v  -  u. 


W  ith  the  above  requirement,  one  can  now  rewrite  Eqs.  (47)  and  (18)  in  terms  of  v 
along;  that  is 


vm 


V  =  — 


m 


(52) 


and 


2P  =  -2—' 


V  m 


(53) 


Equation  (52)  can  be  integrated  to  give 

Vo  V  TT 


Tno 


dm 

m 


(54) 


so  that 


vq  m 


(55) 


and.  defining  the  mass  fraction  ■)  as  m^z  mg 


V  =  Vo 


I  -  ~f 


(56) 


Equation  (56)  is  shown  plotted  in  Fig.  6  as  v/vq  as  a  function  of  y.  For  comparison, 
the  same  velocity  ratio  profile  for  a  constant  1,^  chemical  rocket,  with  u/vq  -  -5  is  shown  as 
a  dashed  curve.  One  can  view  this  result  in  two  ways;  first,  for  the  same  mass  fraction,  the 
hybrid  plume  rocket  can  achieve  significantly  greater  velocities,  and  second,  for  the  same 
terminal  velocity  the  hybrid  plume  rocket  represents  a  substantially  lower  mass  fraction 
and  hence  a  substantiallv  greater  payload  capacity. 


k 


In  this  regard,  it  must  also  be  pointed  out  that  a  constant,  high  I,p  rocket  is  also 
attractive.  However,  the  hybrid  plume  concept  presented  here  has  some  added  advantages 
over  even  this  latter  approach.  This  will  be  shown  next. 

6.2  Rocket  Performance  With  Variable  I,p 

In  order  to  be  desirable  the  hybrid  plume  rocket  must  also  have  sufficient  acceleration. 
■Accordingly,  it  is  important  to  obtain  a  relationship  for  v  as  a  function  of  time.  Such  a 
relationship  will  be  derived  next.  Using  Eqs.  (55)  and  (53)  one  can  write 

2P  dm 

- pdt  =  —  .  (o( ) 

{vomoy  rn^ 


Integrating  this  expression 


(I’omo) 


where  the  negative  sign  has  been  introduced  to  account  for  the  mass  decreasing.  The  result 
after  integrating  is 

-  ..2^  =  i  (59) 

^0  (t’o^o)  ^ 

and  solving  for  m  as  a  function  of  time  gives 


"^(0  --  ,  Jp7-  • 


From  this  expression  one  can  calculate  the  burnup  time  r  which  defines  the  dry  mass  of 


the  rocket  m. 


mo  /  2Pmr 


Equations  (55)  and  (60)  can  now  be  combined  to  obtain  velocity  as  a  function  of  time; 


v(t)  =  Uo 


Como 


and  an  additional  integration  yields  the  distance  as  a  function  of  time: 


x(0  -  Uo< 


I’omo 


m 


'■.w 


Note  from  Eqs.  (62)  and  (63)  that  the  acceleration  is  constant  and  equal  to  'IP/vottiq. 


6.3  Constant  I,p  Rocket  Performance 


By  a  similar  process,  the  conventional  rocket  equation  can  be  shown  to  be 


J’(0  “  i‘o  -  a  In 


1  3t 


(64) 


where  u  is  the  exhaust  velocity  and  3  is  defined  by  the  equation 


3~ 


m 


m/  _  7 


mo  mot  t 


(65) 


and  where  the  mass  flow  rate  m'  is  considered  to  be  a  constant.  The  integration  of  Eq.  (64) 
yields  the  distance  as  a  function  of  time, 


X(t)  -  vot  -  u 


^  -3t  ,  , 

t  -f-  — r —  in(I  -  3t) 


(66) 


and  relating  3  to  power  through  Eq.  (48)  one  ran  use 


3  - 


2P 

u^mo 


(67; 


Equations  (65)  and  (67)  can  be  cf)mbined  to  yield  an  expression  for  the  exhaust 


velocity  or  l,p  in  the  case  where  it  is  a  constant. 


u  = 


2/^ 

mo7 


1/2 


(68) 
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6.4  Illustrative  Examples 


Two  cases  have  been  considered  to  compare  a  variable  /,,,  rocket  to  a  constant  /,p 
one  operating  at  the  same  power  level. 

First,  ccjnsider  a  distance  of  7.8  x  10'  Km.  as  in  a  trip  to  the  planet  .Mars,  and  consider 
twcT  power  levels;  P  .b  MW  and  P  1.0  MW  .  One  inquires  as  to  the  total  trip  time 
for  a  rocket  plus  payload  of  1.122  x  10'^  Kg  delivered  at  the  destination  as  a  function  of  the 
total  fuel  on  board.  I  he  analvsk"  produces  two  times:  t^.  from  F.q.  (66).  which  denotes  the 
trip  time  for  a  constant  I gp  rocket  and  from  Kq.  (62).  which  denotes  that  quantity  for 
a  hybrid  plume  rocket  operating  with  its  exhaust  velocity  optimally  tuned  in  accordance 
with  Eq.  (49).  The  results  are  plotted  in  Fig.  7.  .\<3te  that  the  solution  of  Kq.  (6.3)  gives 
only  one  point  for  a  given  value  of  P  and  hence,  the  two  power  levels  under  consideration 
here  yield  the  two  values  of  t/t  shown  in  the  graph.  Fc^r  on  the  other  hand,  the  solution 
of  Eq.  (66)  yields  a  continuous  curve  simply  because  one  has  flexibility  in  tlie  rn'  and  /,p 
combination  chosen  to  meet  the  power  level  requirement. 

Note  rom  Fig.  7  that  all  the  values  of  tc  are  higlier  than  the  value  of  at  the 
corresponding  power  level.  .Additionally,  decreasing  or  increasing  the  powc-r  level  has  a 
Corresponding  effect  on  trip  time  but  does  not  change  tlie  relative  com[>arison.  Diere  is  a 
minimum  trip  time  and  hence,  an  optimum  fuel  mass  for  the  con -taut  Pp  rocket,  given  a 
[layload  and  a  power  level.  Physically,  this  mean^  that  too  much  fuel  can  significantlv  slow¬ 
down  the  ship  and  too  little  means  a  very  low  acceleration  for  the  trajectory.  The  hvbrid 
plume,  on  the  other  hand,  while  using  more  fuel,  insures  that  the  acceleration  is  constant 
and  trades  off  thrust  for  I ,p  as  the  rocket  velocity  increases. 

The  second  example  is  shown  in  F'ig.  8  for  a  much  longer  trip  and  a  smaller  payload. 
I  he  results  are  qualitatively  the  same:  however,  the  gain  in  performance  of  the  hv  brid 
plume  over  the  constant  I,p  vehicle  is  much  greater  in  this  case. 

These  considerations  show  the  fundamental  advantage  of  variable  /,p  propulsion.  How¬ 
ever.  these  results  must  be  taken  in  the  context  of  rectilinear  motion  in  gravity-free  space. 
The  same  analysis  for  the  gravitational  two-body  problem  of  (jrbit  transfer  is  now  iieing 
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pursued  by  our  research  team.  In  principle,  the  important  cjuantity  to  follow  is  the  vehi¬ 
cle's  stale  vector  (i.e.,  velocity  and  position)  with  respect  to  the  Karth  (if  tlie  fuel  is  loaded 
there).  Under  these  circumstances,  the  relative  velocity  of  the  rocket  may  vary  from  a  low 
value  at  apogee  to  a  high  value  at  perigee  or  gradually  ilecrease  as  the  \'eh.icle  spirals  out. 
The  control  law  for  the  operation  of  the  engine  may  be  different  in  form;  nevertheless, 
variable  !,p  propulsion  will  probably  yield  substantial  advantages  in  this  realm  as  well. 

7.0  Expcriutental  Facility 

The  purpose  of  this  work  is  to  investigate  the  basic  physics  of  the  small  tarulern 
mirror  device.  The  goal  is  to  establish  a  generic  tandem  mirror  plasma  facility  suitable 
for  investigating  the  hybrid  plume  and  other  high  power  x'ariable  /,p  plasma  propulsion 
concepts  1.2.3. 

The  tandem  mirror  design  is  shown  in  Figs.  1  and  3  and  has  the  following  major 
com[X)nents:  the  magnetic  containment  coils,  the  vacuum  vessel,  a  cold  plasma  injector  at 
t  fie  closed  end.  an  rf  antenna  for  ion  cyclotron  resonance  heating  (K’HH)  and  caftability 
for  a  hybrid  plume  exhaust  at  the  open  end.  I  he  plasma  i>  injected  by  a  gun  and  heated  to 
a  few  100  e\  in  the  central  ceil  by  ICfiH.  The  hybrid  plume  tan  be  produced  by  injecting 
a  h\  personic.  coaxial  gas  jet  at  the  exhaust  emi  vvhich  partially  mixi's  'vv  ith  the  hot  [)lasma 
to  prfxiuce  thrust.  The  exhaust  can  be  controlled  fjy  a  (omhunation  of  \ar\mg  the  injector 
current,  the  electric  potential  in  the  anchor,  the  magnetic  field  and  ttas  jet  stagnation 
prc-ssure.  Durirm  tfie  design  pha.se  it  was  decided  to  qiiadrufOe  the  magnetic  field  intensity 
reported  previously  1  tcj  increase  confinement  time.  This  was  clone  v  it h  minor  c  hanges  to 
the  roil  configuration.  The  design  cicf  a  cryc»genic  ccail  svs'tein.  \aruum  vessel  and  structure 
were  kept  as  simple  as  possible.  .Alignment  is  minimal.  Initially,  oiilv  the  ccmtral  cell, 
booster  coil  and  mirror  coils  will  l)e  required.  The  Ioffe  and  saddle  cauls  will  be  added 
later.  Ingenicjus  design  and  fabrication  methods  have  been  devised  to  facilitate  a  drastic 
saving  in  fabrication  costs.  Table  1  lists  the  major  parts  for  the  device,  the  parts  being 
built,  and  the  parts  deferred  for  futiire  upgrade.  The  central  roils  .ind  one  mirror  coil 
have  been  delivered  and  electrical  tests  have  been  performed.  Ilie  wave-form  of  the  tc"-ts 


i>  ''hiiwn  in  Fig.  1.  1  he  iiiirri)r  roil  has  been  energized  up  to  1(,)()()  amp  for  up  to  1  svc  time 
duratit>n  without  cooling.  1  he  true  values  of  inductance  and  resistance  are  alstj  measured. 
The  assembl\’  will  begin  in  the  fall.  1  he  objectives  of  this  e.vperiment  are  listed  in  fable 
11  and  are  also  discussed  beltw. 

laj  IMasina  contaiiuTient.  replenisluTient.  anti  exhaust 

file  iiulk  plasma  will  be  generatetl  by  rf  l)reakdown  of  gaseous  h\  (irogen  in  one  of  tiie 
I  lo>t‘ti  mirror  etui  or  in  the  ttmtral  cell.  1  lie  central  cell  plasma  builtlup  vs  ill  be  obser\'ed 
atid  liie  t'xhaust  tiiroiigh  the  otiier  end  will  be  measured,  funing  of  the  magnetic  field 
toils  anti  electric  plug  potential  prttfile  will  be  performed  tti  obtain  tiie  best  gratle  central 
plasma  and  tiie  most  useful  exhaust.  Means  of  plasma  replenishment  such  as  central  cell 
gas  injectitin.  a  gas  blanket,  pellet  injection,  and  streaming  gas  frtun  the  closed  end  will  be 
investigated.  The  effect  of  these  replenisliing  methods  on  the  exhaust  properties  at  both 
low  and  high  particle  confinement  times,  will  be  the  center  of  interest. 

I  i)  I  Plasma  iieat  ing 

Plug  and  central  cell  plasma  heating  by  ion  cyclotron  resonance  heating  (K'RH)  will 
iu'  Used.  The  efficiency  of  this  energy  injection  method  is  highly  dependent  on  antenna 
design  and  magnetic  field  sha[)ing.  .Accordingly,  a  major  objective  of  this  research  is  to 
optimize  the  antenna  design  for  this  particular  afjplication  and  imestigate  the  operational 
regimes  and  power  levels  which  can  be  achieved.  Other  forms  of  energy  loss,  sudi  as 
radiation  and  charge  exchange  affecting  the  o\erall  power  balance,  will  be  measured  as 
well. 

(c)  Stability  and  cfintrol 

One  of  the  major  objectives  of  this  experiment  is  to  generate  stable  and  well-behaved 
plasmas  over  as  large  a  set  of  operating  coiiditiirns  as  possible  for  the  compact  device.  A 
great  deal  of  existing  experimental  knctwledge  and  raw  data,  derived  from  similar  fusion 
experiments  will  be  utilized  here  to  achieve  the  best  plasma  operation  under  \arying  con¬ 
ditions.  I  bis  will  become  useful  later,  in  future  experiments  with  variable  /,p  exhausts, 
f dearly  then,  stability  and  control  of  the  plasma  will  be  an  important  area  of  investigation. 


Ml  )rei  the  \ery  fact  that  a  preferential  flow  direction  exists  in  this  case  will  mean  a 
certain  asymmetry  in  the  configuration  of  the  machine  iri  the  longitudinal  direction.  I'lie 
effects  of  such  asymmetries  on  plasma  stability  will  be  investigated. 


1  he  bulk  of  the  experimental  ilata  will  be  integrated  in  the  form  of  an  "operational 
eiuelope"  map  for  this  device.  Such  a  map  will  continue  ti)  be  ref..,ed  as  more  experience 
i'  gained  and  will  form  the  basis  for  the  next  phase,  where  a  wide  variety  of  high  thrust, 
xariable  rocket  ex[)CTiments  will  fie  invest igi.: cd. 

id)  Instrumentation 

(X)ntrol  system  will  be  developed  to  achieve  high  reproductivity  of  the  plasma  and 
exhaust  properties  with  minimum  monitoring  for  reliable  operation. 
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TABLE  I 


Major  Parts 


Parts  Being  Built  Parts  for  Upgrade 


Central  Cell  Coils 

X'acuum  Chamber 

Booster  Coils 

Saddle  Coils 

Ioffe  Coils 

RF  Transmitter 
and  Antenna 


Central  Cell  Coils 
\'acunm  Chamber 
Booster  Coils 


Saddle  Coils 
Ioffe  Coils 


RF  Transmitter 
and  Antenna 


TABLE  I  I 


Experimental  Objectives 

1.  Plasma  containment. 

RF  brecikdown  of  hydrogen  gas. 

Central  plasma  buildup. 

Scanning  of  magnetic  field  and  electric  plugging  field  for  best  grade  central 
plasma. 

2.  Plasma  Heating. 

Plug  and  central  cell  plasma  heating  by  ion  cyclotron  resonance  heating 
method. 

Optimize  heating  efficiency. 

3.  Stability  Control. 

Generate  stable  and  well-behaved  plasmas  over  a  large  set  of  operating 
conditions. 

4.  Instrumentation 
Minimize  instrumentation. 

5.  Thrust  and  I,p  study. 


z-axis  (cm) 


z-axis  (can) 


^51  CCNTOUR  PLOT 


P51  CONTOUR 


z-axis  (cm) 


z-axis  icm; 


Figure  2 

\Ia^netlc  field  lines  with  neutral  i;as  jet  near  wall  (svnunetry  line  at  the  r..p  ileach  pi  jt  i  Ni-utral 
cas  IS  "ionized"  at  ■10cm  down  duct,  and  momentum  pushes  field  lines  plasma  iluwn  duet  The  flow 
in  these  figures  is  not  self  consistent,  but  represent  the  ideal  fluid  flow  possible  ilai  Initial  field 
configuration,  (lb)  Island  (or  smoke  ring")  begins  to  form  ilc)  Island  breaks  awav  ildi  Island 
continues  to  move  down  d  iet,  disconnected  from  nozzle 


(a) 

5  Volts 
220  amp 
1  sec 


(b) 

2  Volts 
82  amp 
1  sec 


Inductance  3  mh 

Resistance  (room  temp)  30  m  ohm 

Figure  4:  Typical  voltage  (upper  curve)  and  current  (lower  curve) 
traces  for  energizing  the  central  cell  coil  at 
(a)  220  amp  and  (b)  82  amp. 
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Total  payload  is 
1.122x104  Kg. 
Assume  total  trie 


Assume  total  trip 
distance;  7.8x1  O^Km 
Tc  =  constant  Up  time 
Th  =  variable  Up  time 
P  =  power  level 


Tc(mJatP=  .5  MW 


Th  at  .5  MW. 
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